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In this letter we study both ground state properties and the superfluid transition temperature of 
a spin-1/2 Fermi gas across a Feshbach resonance with a synthetic spin-orbit coupling, using mean- 
field theory and exact solution of two-body problem. We show that a strong spin-orbit coupling can 
significantly enhance the pairing gap for l/(fcFas) < due to increased density-of-state. Strong spin- 
orbit coupling also significantly enhances the superfiuid transition temperature when l/{kpas) < 0, 
while suppresses it slightly when l/(A;Fas) ^ 0. The universal interaction energy and pair size at 
resonance are also discussed. 



During the last few years, studies of ultracold Fermi 
gases across a Feshbach resonance (FR) have brought a 
lot of excitements to physics [1]. On a separate devel- 
opment, recent experimental breakthrough on synthetic 
gauge field has open up a lot of new opportunities to cold 
atom physics [2, ^31. One application of this technique is 
to engineer an effective spin-orbit coupling (SOC) in cold 
atom system [4]. Very recently, a pioneer experiment in 
NIST has already achieved a restricted class of spin-orbit 
coupled BEG of *^Rb atoms [3]. Theoretically, even the 
mean-field study of boson condensate with SOC has re- 
vealed many interesting physics [SHZ]- For fermions, a 
concrete scheme has also been proposed for generating 
SOC in ^'^K atom in the regime where a magnetic FR 
is available [S], and the experiment of implementing this 
proposal is now going on in the laboratory. However, a 
theoretical study of spin-orbit coupled Fermi gases across 
a FR is still lacking. 

In the absence of SOC, a Fermi gas across a FR pos- 
sesses three key physical properties: i) across a FR the 
system undergoes a crossover from a BCS type fermion 
superfluid to a BEC of molecules; ii) at the FR, it is a 
strongly interacting system and exhibits many universal 
behaviors; iii) nearby a FR, the transition temperature 
of fermion superfluid T^/T-p is the highest one among all 
fermion superfluids (or superconductors). The question 
is that how these three properties evolve in the pres- 
ence of SOC. For i), since now the pair wave-function 
has more complicated structure with both singlet and 
triplet components, and exhibits p-wave character in the 
helicity bases, one needs to investigate whether it is still 
a crossover or there is a phase transition in between. 
Even if it is still a crossover, how SOC affects it. For 
ii), since the strength of SOC introduces another length 
scale X/kp, the universal constants at resonance now be- 
come universal functions of A/fcp, and we want to un- 
derstand the behaviors of these functions. And for iii), 
the question is whether Tc/Tp will be enhanced or sup- 
pressed by SOC. (Here the units kp, Ep and Tp are the 
Fermi momentum, the Fermi energy and the Fermi tem- 
perature for non-interacting system without SOC. as is 
the s-wave scattering length.) 

In this letter we address these issues using both mean- 



field (MF) theory and exact solution of two-body (TB) 
problem, and the main results are summarized as follows: 

(Al) The system is gapped for all values of Ug. The 
pair wave-function obtained from MF theory has the 
same symmetry property as the wave-function of TB 
bound state, and they coincide with each other for 
l/(fcFas) S> 1. These two evidences support a crossover 
picture instead of a phase transition. 

(A2) The order parameter A always increases as the 
strength of SOC X/kp increases. For l/(fcFas) ^ 0, the 
increasing becomes profound when X/kp is large enough 
that the density-of-state (DOS) at Fermi surface is signif- 
icantly enhanced. While for l/(fcFas) ^ the increasing 
of A is always less significant when the chemical potential 
drops below the single particle energy minimum. 

(B) At resonance, the interaction energy Ei^t/Ep and 
the pair size kpl as functions of X/kp have very different 
behaviors for X/kp <C 1 or ^ 1. 

(CI) For l/(fcFas) < 0, is enhanced by SOC due 
to two effects. One is that the increased DOS enhances 
Tbcs from MF theory; and the other is because stable 
molecules with finite binding energy now also exist in 
this regime. 

(C2) For l/(fcFas) > 0, is given by BEC temper- 
ature of molecules Tbec, and is slightly suppressed by 
SOC because the effective mass of molecules is increased. 

(C3) At resonance Ug — ±oo, Tc will finally saturate 
to 0.193Tf when X/kp is large enough, which is higher 
than the transition temperature without SOC. 

Model: We consider an isotropic in-plane SOC. The 
single particle Hamiltonian is given by Hq — Tp^ /{2m) + 
Ap_L ■Tj./m, where p_L = {Px,Py) and cr_L = icrx,cry) (set 
A > and h = 1). The generalization to anisotropic and 
more complicated SOC is quite straightforward. In the 
second quantized form, Hq ^ I]p[ep(Cp.,.Cp^ -I- cJ^^Cp^) -I- 

Ap_L(e"'^pc^.,.Cp^ + e^'^pc^^Cp^)], where Cp = pV(2™), 
P± = IP-l| and ipp ~ aj:g{px + ipy)- The single particle 
Hamiltonian can be diagonalized in the helicity bases as 

"^0 = Ep[Cp+'ip,+^P,-(-+Cp-^p,-^P,-]> with^pi = ep± 
Apj_/m, where helicity ± means that the in-plane spin is 
parallel or anti-parallel to the in-plane momentum. The 
fermion operators in the helicity bases are related to the 
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FIG. 1: (a,b) Order parameter A/Ef and chemical potential 
fi/Ep (measured from single particle minimum) as a function 
of l/{kpas) for different X/kp; (c) A/Aq as a function of A/fcp 
for three different values of l/(fcpas), where Ao is the paring 
gap without SOC. (d) DOS at Fermi energy A/'(/i) compared 
to DOS without SOC (A/'o(i?p)) as a function of A/fcp. 



fermion operators in the original spin bases via /ip,+ = 
(cpt + e-*'^pcp;)/V2 and = (e^'^pcp^ - CpJ/\/2. 

When the effective range vq of inter-atomic poten- 
tial is much smaller than all the other length scales 
in the problem, i.e. fcp^'o ^ 1 and Xtq <C 1, as 
in the conventional crossover theory, we use a zero- 
range potential to describe the intereaction between 
atoms. The interaction can be written as Hmt — 

(g/y) Epp'q Cq/2 + ptCq/2-piCq/2-p'iCq/2 + p't, whcrC g iS 

related to via 1/g = mj (47ras) — l/(2ek), and V is 
the system volume. 

Mean-field Theory: For MF discussion, we only focus 
on q = channel of H-mt- Introducing the order param- 
eter A = (g/V) X]p(c-p^Cp^), one can obtain the mean- 
field interaction 7^?^^ = AJ2p{clfclp^ + 'ti.c.)-\A\^V/g. 
Transforming it into the helicity bases, it becomes 
Hf^^ - -(A/2)Ep(e-^'^-/ii,,+/i^_p.+ + e^^-/zi.,_/i^_p,_)- 
|ApV^/5. One can see that pairing only exits between 
atoms with same helicity, and the pairing of helicity ± 
has Px T 'iPy symmetry. 

Hence, the MF Hamiltonian is given by Hmf — "Ho + 



Hf^^ — fiN. It is very to solve Hmf in the helicity bases, 
which gives 
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where ep± = y (Cp± — /i)^ -t- is the energy of quasi- 
particles and fp^± = tanh[ep±/(2fcBr)]. 



Two-body Problem: The TB problem in the presence 
of SO coupling has been solved in Ref. f9^ for the case 
of molecular center-of-mass momentum q = 0. It was 
found that the TB bound state appears even at the BCS 
side of resonance with < 0, because of the increase of 
low-energy DOS [S] . Here we solve the two-body problem 
for finite q, which is very useful for later discussions. In 
general, the TB wave-function can be assumed as |^)q = 

EL[V'n(k)c|j/2+ktCq/2- 
V'n(k)c|j/2 + ktCq/2-kt - 

means the summation is over half of momentum 
space. The Schrodingcr equation {Hq -\- Hint)|^')q = 
^'q|^)q leads to a self-consistency equation as [1^ 



■q/2H 
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^u(k)Cq/2 + kiCq/2 
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where fk.q = ^^q - eq/2-hk " eq/2-k, and (^^kq = ipk - (Pq- 

For q = 0, Eq. ^ recovers the results in Ref. [5], and 

for any there is always a bound state solution Eq < 
— A^/m. We obtain an analytical equation for Eq 



2^/{-Eo)m- Xln 



y/{~Eo)m + X 
^{~Ea)m-\' 



(4) 



With i?o, one can then use the Schrodinger equation to 
determine the bound state wave-function. 

Results and Discussions: With the MF theory and TB 
solution presented above, we are ready to address the 
questions posted at the beginning. 

(A) Crossover: We solve the MF equation for T — {). 
The value of order parameter A / E-p and chemical poten- 
tial /i/i?F as a function of l/(fcFas) is shown in Fig. [ija- 
b), for various A/fcp. Not surprisingly, their behaviors are 
not qualitatively different from A = 0. Nevertheless, it is 
worth to point out that for conventional p-wave pairing, 
even though the order parameter is non-zero, the pair- 
ing gap will still close when the Fermi surface touches 
Pj^ = point. However, the paring form factor here is in 
fact [px ± ipy ) Ipjl instead of conventional px ± ipy , which 
ensures that the superfluid is always gapped. 
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FIG. 2: (a) rj — ($|»I')q=o as a function of l/(A;Fas) for differ- 
ent A/fcp. (b) The ratio of triplet component to singlet com- 
ponent Zt for TB wave-function and MF pair wave-function, 
for two different A/fcp. 
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From MF theory, we can obtain the BCS wave-function 



|BCS) cx exp 
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where Uk.i = e^'^Vk 
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function as |<i>) = 

0tt(k)cLc;Lkt ^ 



and itk^i = 
Then we can define a pair wave- 
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'?^u(k)4iC^-ki]|0)' where 



^/c ( 




+ 


Wk - 


) 




"k,+ 









•^kH 

Mk,- 



^k- 
Wk,- 



(5) 
(6) 



and C is the normalization factor. The symmetry prop- 
erties of pair wave- function | $) agree with that of zero- 
momentum molecular wave-function |^')q=o discussed in 
Ref. [9j. We then compute their overlap rj = ('I>|^')q=o 
as a function of l/(fcFas), as shown in Fig. [2]ja). The 
overlap approaches unity rapidly when l/(fcFas) ^ 0. 
In Fig. [2][b) we also plot the ratio of triplet to singlet 
component for both |$) and |^')q=o- It shows that the 
TB wave-function always has a larger triplet component. 
Nevertheless, they converge together quickly. 

In Fig. [ijc) we plot A as a function of A/fcp, from 
which one can see that there is a characteristic value 
roughly located at A/fcp ~ 0.5. Below this value the 
change of A with A/fcp is small, while above this value 
the increasing of A becomes very significant. In Fig[l|d) 
we show the DOS at Fermi energy N{^) compared to the 
DOS without SOC (A/'o(i?F))- Their ratio remains nearly 
unity until reaching A/fcp w 0.5, and then it increases 
rapidly. It is because for low density or strong SOC, 
the Fermi energy drops below Dirac point at pj^ — 
with fcp < (37r/4)-'^/^A, and only the lower helicity minus 
branch will be occupied. In this case its DOS N{^) — 
mA/(27r) is a constant independent of ^, while without 
SOC, N{£^) ^ -yj, therefore the DOS is always increased 

by SOC nn]. 

By comparing Fig. [iJc) and (d) one can draw the con- 
clusion that the increasing of A is due to the increasing 
of DOS. From Fig. [3a,b), we also notice that when /x 
decreases below the single particle energy minimum, the 
DOS effect is no longer important. Thus, the influence 
of SOC on A becomes very weak. 

(B) Universality: As we all know very well now, at 
resonance when — >■ ±oo, the interaction energy per 
particle Eiat/Ep will not diverge, instead, it saturates to 
a universal value of the order of unity. Now, this universal 
value becomes a function of A/fcp. Within MF theory, 
we define the interaction energy as i?int = ((BCSlHo + 
?iint|BCS) — Eq)/N, where Eq is the total energy of a 
non- interacting system. In Fig. [sja), we plot Ei^t/Ep 
as a function of A/fcp at resonance. Its behavior is very 



different in the regime of small and large A/fcp. For A <C 
fcp, we have Ei^it/Ep w —0.24 -I- o{X/kp); while for A >• 
fcp, we find Eint/Ep ~ — 0.44(A/fcp)^. It is because from 
Eq. ^ one can find out that at resonance, TB bound 
state energy Eq = — 2.88A^/(2m), and the binding energy 
is given by — A^/rn — Eq = 0.88A^/(2m), which is twice 
of — E'int in the limit of strong SOC. 

Another notable feature of unitary regime is that the 
size of Copper pairs fcp^ is also of the order of unity. Here 
we can compute the anisotropic pair size from 




where a — x,y and z. When A/fcp ^ 0, Ix = ly < Iz 
which means the Cooper pairs are elongated, as shown 
in Fig. [Sjjb). Similarly, we find for small A/fcp, fcp^ ~ 
1 + o(A/fcp), while for large A/fcp, fcpZ cx fcp/A. Their be- 
haviors at large A/fcp shows that the system still behaves 
like weakly interacting molecular BEC. 

(C) Superfuid Transition Temperature: With 
MF theory, one can calculate the BCS temperature 
Tbcs/Tf as shown in Fig. |4j It increases as A/fcp in- 
creases, for the same reason of DOS effect. 

SOC also affects the effective mass of molecules. Sub- 
stituting the molecular dispersion Eq — Eq + q'^ / {2mi,) + 
ql/{4m) into Eq. and expanding Eq. ([sj to the order 
of q^, we obtain an equation satisfied by as 
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where £k,o = Eq — 2ek, and Eq as a function of l/(Aas) 
can be obtained from Eq. Q. Solving Eq. ^ one can 
find a relation between rrih and Eq as [lOJ 
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Since Eq < —}? jm < 0, mb is a monotonically decreas- 
ing as l/(Aas), as shown in the inset of Fig. |4j Since the 





FIG. 3: (a) Eint/Ep as a function of A/fcp. The dashed line 
is a fit of — 0.44(A/fcF)'^; (b) Size of Cooper pair in x-y plane 
Ixkp and along z direction Izky as functions of A/fep. Both 
are plotted at resonance for Os = txj. 
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FIG. 4: Superfluid transition temperature Tc/Tf from 
MF theory (Tbcs dotted lines), from BEG temperature of 
molecules (Tbec dashed lines), and from an interpolation 
scheme including the contributions from non-condensed pairs 
(solid lines). The symbol with error bar is the Monte Garlo 
results without SOC [11]. Inset: The in-plane effective mass 
of TB molecule rrib/m as a function of l/(Aas). 



bound state always exists for any as 0, we can discuss 
the BEC temperature of molecules at both sides of res- 
onance, which is given by Tbec/Jf = 0.218(2m/mb)^^^. 
When l/(Aas) ^ — oo, 771b — 4772, which gives Tg^Q — 
O.ISTTf- For a given negative Og, TTib decreases as A in- 
creases, and thus Tbec increases. When l/(Aas) — >■ +oo, 
nih = 2777, which gives Tbec = 0.218Tf. For a given 
positive flg, 7r7b increases as A increases, and thus Tbec 
decreases. At resonance — oo, mi,/m = 2 AO is a uni- 
versal value, and one obtains Tbec — O.lOSTp. When 
A/fcp is large enough that the molecules become tightly 
bound, the actual Tc should be very close to TbeCj which 
is higher than Tc = O.lSTp without SOC ^llj. 

A controllable calculation of superfluid transition tem- 
perature in the entire crossover regime is a difficult task 
even without SOC. A widely used approximation scheme 
is the NSR method [H], in which Tc is determined by 
Thouless criterion and a modified number equation 



771 

47ras 



AV 



/p 



(^p,+ -m) (Cp,--Ai) Cp 



77 = 77flu 



-E 



V 



fp,+ fp. 



(8) 
(9) 



where /p_± = tanh[(^p_± — /z)/(2fcBTc)]. The number 
equation contains the contributions from free fermions 
and non-condensed bosonic pairs n^uc- In the NSR ap- 
proach, 7ifluc can be obtained from diagrammatic calcu- 
lations. However, such a calculation becomes much in- 
volved in the presence of SOC, and we leave it for future 
investigations. Here, as a rough estimation, we interpo- 
late Tc between Tbcs and Tbec by making the approxi- 
mation 77fluc = i^pl/[e((fl/(2m,)+p^/(4m))/(fcBT<.) _ 

In fact, such an approximation is quite reasonable in 



crossover regime for A/fcp ^ 1, since as one can see from 
Fig. [sjjb), the size of pairs is already smaller than inter- 
particle distance. The interpolation results are shown as 
the solid line in Fig|4] from which one can see Tc is signifi- 
cantly enhanced for l/{k-pas) < 0; while for l/(A:Fas) > 0, 
the suppression is insignificant. 

As an initial effort to understand this rich system, this 
work points out some basic features as summarized at the 
beginning with simple techniques, and leave more accu- 
rate studies with more advanced techniques for future in- 
vestigations. Our predications can be verified experimen- 
tally once such a system is realized. These studies are also 
first step toward interesting topological phases in this sys- 
tem with population imbalanced in two-dimension. 
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Appendix: In this supplementary material, we present some details of calculating the density of state and solving the 
two-body Schrodinger equation. 

Density of State 

From the single-particle dispersion ^p,± = ep ± Xp±/m, we obtain density of state for each helicity branch as 



m 

m 
27r2 



— T ttA . A 
' 2mt h A arctan , 

/- — - ttA A 

' im4 H h A arctan 

2 V^rn^ 



2m'' 



(10) 
(11) 



where is the step function. Comparing to the case without SOC, total DoS Af{^) = A/"- (^) -|-A/ + (^) is always enhanced. 
For — ^ < ^ < 0, only Af-{$,) contributes to the total DoS which is a constant independent of energy. 
In a ideal Fermi gas with particle density n, chemical potential at zero temperature is determined by 



(12) 



For n < A''/(47r^), only the helicity minus branch is occupied, and we find /u = |^Sf — Hence in this case, DoS at 
Fermi surface is increased linearly as a function of SOC strength A, 

^(-"^ = > (3^)'^'^^) (13) 

where M)(£^f) = mk^/'K^ is the DoS at Fermi surface without SOC. 



Solution of Tow-body Problem 

The two-body wave-function with a center-mass momentum q can be written as 

= E' [V'n(k)c|+k,4_^^ + V^t(k)4^^^4_^^ + Vn(k)4^^^4_,, + Vu(k)4^^,4_,J |0), 

k 

where denotes a summation with kz > 0. Schrodinger equation for two-particles interacting via a contact potential 

(7^0+ Hint) |*)q = Sq|*)q (14) 

can be written explicitly as 

[^n(k') - ^,t(k')] + ^ [(t - fe=^) + ^(f - ky)] V'n(k) + ^ [(f + k.) - i{\ + fe.)] Vu(k), 

k' 

; [(f - fe) - «(f - fc.)]^n(k) + ^ [(f + fe) - i(f + fc.)]^u(k), 
■ [(f + A;,) + i(f + fc,)] Vn(k) + ^ [(f - k.) + i(f - fc,)] Vit(k), 



fk,qVn(k) = 



fk,cVit(k) = K(k') - Vn(k')] + ^ [(f + '^x) + »(f + fc.)] V'tt(k) + ^ [(t - fc-) - »(f - fey)] V'u(k), 



fk,qV'n(k) 

fk,qV'u(k) 
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with fk.q = - eg+k - <;s_k- 



Introducing (,6s(k) = ^ ['(,/;j^^(k) — -(/;^^(k)], '!/'a(k) — [V't-I. (k) + V'.|,t(k)], which are corresponding to the wave-function 
components | ti — it) y-ucl I ti + it) respectively, one can find a self-consistency equation for -i/is 



m'^Sk^q m2£'k,q(m2f^ q 
Hence, the energy eigenvalue Eq is determined by 
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(16) 



whore t^qk = '/'q ^ "i^k is the in-plane angle between q and k. With obtained Eq, the TB wave-function can also be 
determined through Schrodinger equation. 
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Bound State with Zero Center-of-Mass Momentum 



For the eigen-state with center-of-mass momentum q = 0, the threshold energy of the scattering state is — A^/m, and the 
energy of molecular bound state is determined by 
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Introducing the molecular binding energy Eh = — A^/m — Eo, above equation can be re- written as 

1 _ VmEb + A2 _ 1 VmEb + A^ + A 
Afls ^ A 2 y^mE^, + A2 - A ' 



(17) 



(18) 



One can see that for any value of as there is always a bound state with i5b > 0. In the limit of Atts — )■ 0^, |iJo| 3> A^/m, 
we obtain 



Eu 



(19) 



which is same as the case without SOC. In the limit of Aas — s> , |i5o| — A^/m, we obtain 

4A^ _ 2 

= — e 

which shows a exponential dependence on (Aas)~^. At unitarity, as — >■ oo, we find 

,A2 



£b = 0.439- 



The bound state wave-function is given by 

V't-t(k) = -V'it(k) 
V'n(k) = -'0u(k) 
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-Eo — Ck,+ Eo — i^k,- 



(20) 
(21) 

(22) 
(23) 



where C' is the normalization coefficient. The fact ?/'a(k) = implies that for the bound state with q = the triplet 
component of | + it) vanishes. 



Effective Mass of Molecule 



For the bound state with a finite center-of-mass momentum, the eigen-energy can be written as Eq — Eo + q^/{2mh) + 
q^/{Am) if q is small enough. Substituting this dispersion in Eq. ( 16 1 and expanding to the order of q\, we find 



2m ^ £k,o + 



rA'fei 



16A sin ^p\^ 



rA2fc2 



y — 



where fk,o ~ Eo — k^ jm. The integral on the l.h.s. above can be computed straightforwardly as 



-y 



V 



rA2/cr 



-y 



2V 



1 



+ 



1 



m 

8^ JO 
167r JO 



(fk,o-^Afcx)2 ' (fk,o + ^Afcx)2 



dk_ 
dk 



dkz 



+ 



ki 



(mEo - fc2^ - fc2 - 2Afcx)2 (m-Eo-fci-iti-f 2Afcx)2 

k± k± 



{k\ + 2Xk± - mEo)^/2 {kl - 2Afcx - mEo)'-^^^ 



^Jm\Eo\ 
8^m|£o| - A2' 



(24) 



(25) 



7 



and the integral on the r.h.s of Eq. ( 24 1 can be computed as 



-T 



-T 

k 

2 \ poo 

m A ' 



A sin^ ((9k 



V ^ 2k^ 

k 



1 



1 



167r2 



m^A 
' 32-K 



dfci 



(fk.o - ^Afcx)^ (fk.o + ^Afcx)2 AfcxVfk.o-I^Afcx fk.o + J^Afc^ fk.o 
dfc. 



(mSo -kl-kl- 2Afcx)2 (m£o - fc^^ - fef + 2Afcx)2 
1 1 



vfc, V 



+ 



Afcx \mEo — k\ — k'i — 2\k± mEo — k'j^ — k'i + 2\k±_ mEo — k'] — k^, 



dkA 



{k\ + 2Afcx - mEo)'^/^ (k\ - 2Afcx - mEo)^/^ 

2/1 1 



16tt ^/rn\Ea\ 



In 



Afcx V ^fe2 2Afcx - mEo ^ \/k'i - 2\kj_ - mEo \/k\ - mEo 
m\Eo\ A^ 



\Eo\-X^ m\Eo\-\^ 



Thus the in-plane molecular effective mass mb is given by 



2m _ ^ 1 
mb 2 

= 1-1 



ml^ol - A^ ml^ol - A^ A^ 
mj_Eoj m\Eo\ m\Eo\ 

mE]j J mBb A^ 

m£;b + A2 ^ mEi, + A^ ^ m£b + A^ 



One can see that mb is always larger than 2m for non-zero A. In the limit of Aos — )■ 0^, i?b 2> A^/m, we find 

mb = 2771. 

In the limit of Aos — >■ 0~, i5b ^ A^/m, we find 



At unitarity, i5b ~ 0.439A^/m, we obtain 



mb = 4m. 
mb = 2.40m. 



(26) 



(27) 



